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Abstract 

Superconductivity (S) and ferromagnetism (F) are probed through transport and magnetization 
measurements in nanometer scale HoNis-NbN (F-S) bilayers and HoNis-NbN-HoNis (F-S-F) trilay- 
ers. The choice of materials has been made on the basis of their comparable ordering temperatures 
and strong magnetic anisotropy in H0M5. We observe the normal state reentrant behavior in 
resistance vs. temperature plots of the F-S-F structures just below the superconducting transition 
in the limited range of HoNis layer thickness djjN (20 nm < duN < 80 nm) when d^bN is fixed 
at ~ 10 nm. The reentrance is quenched by increasing the out-of-plane (Hj_) magnetic field and 
transport current where as in-plane (Hm) field of < 1500 Oe has no effect on the reentrance. The 
thermally activated flux flow characteristics of the S, F-S and F-S-F layers reveal a transition from 
collective pinning to single vortex pinning as we place F layers on both sides of the S film. The 
origin of the reentrant behavior seen here in the range of 0.74 < Tcurie/^C < 0.92 is attribute 
to a delicate balance between the magnetic exchange energy and the condensation energy in the 
inter facial regions of the trilayer. 

PACS numbers: 74.45. +c, 74.62.-c, 74.25.Wx 
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INTRODUCTION 



The antagonistic order parameters of superconductor (S) and ferromagnet (F) lead to 
several fascinating effects in the transport and magnetic properties of thin film S-F het- 



[i. 



erostructures HI, |2] . It has been noticed that there is a large suppression of the supercon- 
ductivity at the S-F interface due to the strong pair breaking effect of the ferromagnet via 
spin-flip scattering and/or spin rotation. The ferromagnetic layer is also affected by the 
presence of the superconductor as the Cooper pairs entering the F region acquire a center 
of mass momentum due to the exchange field of the F. This adds an oscillating term to the 
Cooper pair wave function inside the ferromagnetic region {^J . The experimental studies on 
S-F bilayers ^-7|, trilayers jsiOl multilayers 11- 14 1 and S-F-S junctions 15-20] show ef- 
fects such as the direct and inverse proximity effect, reentrant superconductivity and critical 
temperature oscillations. 

In addition, the phenomenon of long range triplet pairing is possible when the coupled 
ferromagnetic layers have inhomogeneities in magnetization, which can be due to domain 



walls, spiral magnetism and spin scattering at the interface [19|-|2l|. In the case of a F-S- 
F trilayer, for certain angles between the magnetization vectors of the F layers, a triplet 
pairing can be induced in the F layers 23|, |5lj. As a result, the Cooper pair can survive in 
ferromagnetic region up to length scale of normal metal coherence length (£jv)- 

In most of the S-F hybrids studied till date, the F order sets in at a temperature (T Cu rie) 
higher than the transition temperature (Tc) of the superconductor and thus the nucleation 
of the superconducting state is subjected to a robust exchange field of the ferromagnet. 
An equally important and perhaps much more illuminating option is to have a ferromagnet 
whose Tcurie is lower than the Tc of the superconductor. If the ferromagnet is sufficiently 
thin, one may see suppression of Tcurie due to induced superconducting order in the F. To 
the best of our knowledge, there is no direct evidence of this, while the converse has been 
well demonstrated. 

The heterostructures of NbN and H0M5 offer an ideal system to study suppression of 
F order due to superconductivity because NbN becomes superconducting at ~16 K and 
the bulk ferromagnetic ordering temperature (Tcurie) of H0M5 ~5 K. In our previous ar- 
ticle on NbN/HoNis bilayers 24J, we have reported a spectacular flux flow induced peak 
effect in the magnetoresistance R(H) of the system at temperatures below which HoNi 5 is 
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ferromagnetic. This observation was attributed to a spin reorientation transition induced 
by the in-plane magnetic field. These bilayer samples, however, do not show any direct 
signatures of the competition between the F and S order parameters. In contrast, a trilayer 
film of HoNi 5 /NbN/HoNi 5 displays dramatic effect of the antagonism in R(T) and M(T) 
measurements as the system is cooled through the superconducting and magnetic transition 
temperatures. The most striking signature of the competition is a reentrant superconducting 
transition in the R(T) data of the trilayers when the NbN films are made sufficiently thin 
such that the perturbation from the F layer is strong enough to suppress superconductivity 
over a limited range of temperature below the T^ of NbN. We have established a critical 
range of Tcurie/Tc, where Tcurie and T c respectively are magnetic and superconducting 
transition temperature of the trilayer in which the reentrance is seen. The reentrant behav- 
ior is also seen prominently in the M(T) data. We draw a comparison between these results 
and a recent prediction of reentrant superconductivity (RES) in F-S layers where F layer 
has inhomogeneous magnetization 25[. Since H0M5 is a highly anisotropic ferromagnet 
and the films used here are polycrystalline with evidence of inhomogeneous magnetization, 
it may allow a strong superconducting proximity effect at lower temperatures. However, 
we can not rule out the possibility that Ho 4f moments at both interfaces perturb the su- 
perconductivity of the thin NbN layer, leading to reentrance to the normal state. To the 
best of our knowledge we are presenting first observations of reentrant superconductivity 
in S-F heterostructures of Tcurie < Tc, while it has been observed prominently in bulk 
superconductors such as ErRh 4 B 4 , HoMoeSs and HoNi 2 B 2 C [261-130 1. 



EXPERIMENTAL 



Thin films of NbN, H0M5 and their heterostructures were deposited on (100) cut MgO 
substrates using pulsed laser deposition technique. For HoNi 5 , a stoichiometric polycrys- 
talline target was ablated in 0.1 mbar neon environment and then a NbN film was deposited 
in 0.1 mbar N 2 pressure on top of the HoNi 5 layer. For the deposition of the top HoNi 5 
layer, the nitrogen was flushed out and the growth was carried out at 0.1 mbar of neon. The 
growth temperature for all three layers was 200° C. The thickness dffN of the bottom and 
top H0M5 layers was kept the same but varied from 20 to 80 nm in different samples and 
the NbN thickness (djy&jv) was changed from 10 to 30 nm for given <1hn- Further details 
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of growth condition can be found in previous reports [3l|, |32|]. For the measurements of 
electrical transport, the samples were patterned into 15 /im wide lines by using tungsten 
shadow mask and Ar + ion milling technique. A schematic diagram of sample geometry is 
shown in the upper inset of Fig. 1(a). 

RESULTS AND DISCUSSION 

We first present the results of x-ray reflectivity measurements on HoNi 5 -NbN bilayers of 
d-HN — 50 nm to quantify the extent of interface roughness. These data are shown in Fig. 
1(a). A best fitting of a genetic algorithm (solid line) yields a roughness of 0.7 and 1.3 nm for 
NbN/HoNis bilayer with 30 and 10 nm thick NbN respectively. The superconducting transi- 
tions as seen in resistivity measurements on HoNis-NbN-HoNis trilayers of NbN thickness 10 
and 30 nm are displayed in the lower inset of Fig. 1(a). The trilayer with 30 nm NbN shows 
a sharp transition with critical temperature (Tc) of ~11 K. On reducing the thickness to 
~10 nm, the onset of Tc shifts down to ~ 6 K. This drop, however, is truncated by a valley 
and then upturn of R(T) on reducing the temperature below ~ 5.5 K. The resistance goes 
through a peak at w 5 K followed by a sharp drop at lower temperatures. This behavior 
is similar to the reentrant phenomenon seen in bulk samples of HoNi 2 B 2 C, ErRh 4 B 4 and 



HoMo fi S 



6^8 
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3|-[30|]. A magnified view of the reentrant nature of the R(T) of the trilayer with 
10 nm thick NbN is shown in panel (b) of Fig. 1. The R(T) is separated into three regions: 
a normal metallic behavior for T > 6 K, reentrant superconductivity (red bar) between 5 K 
< T < 6 K and zero resistance state for T < 5 K. No such reentrant behavior was seen in 
the trilayers with 30 nm thick NbN. In the inset of Fig. 1(a), we also display the R(T) data 
for a bare 10 nm thick NbN film. From a comparison of curve A and B, it is clear that the 
Tc is suppressed drastically in the trilayer for a fixed NbN layer thickness. 

At this stage it is important to recall the magnetic properties of HoNi 5 , which shows bulk 
ferromagnetic ordering at ~5 K. The magnetism in this compound is highly anisotropic 
due to 4f electrons of Holmium. The easy axis of magnetization lies in the ab-plane of 
the hexagonal unit cell, but the magnetization for field along the c-axis does not saturate 
at field as high as 15 T 33] . The HoNi 5 thin films grown on (100) MgO investigated by 



our group 



31] are polycrystalline in nature with grain size in the range of ~ 10 nm. To 



establish magnetism and superconductivity in HoNi 5 /NbN/HoNi 5 trilayer with 10 nm thick 
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NbN, we have measured temperature dependent magnetization M(T) after zero field cooling 
the sample and then applying 100 Oe field along the out-of-plane direction. The trilayer 
used in magnetization measurements was grown along with the sample used for transport 
measurements in the same run. The M(T) shown in Fig. 1(c) rises sharply below 5 K, 
goes through a peak and then drop to a negative value, indicating strong diamagnetism. 
A magnified view of M(T) in the temperature range of 3 to 7 K is shown in the inset of 
the figure along with the M(T) of a bare HoNis film on MgO [31]. The Tcurie extracted 
from the behavior M in the crtical regime ~ Mq(1 — T/Tcurie) 13 is ~ 5.5 K and ft ~ 0.52. 
As we see in Fig. 1(b) the onset of superconductivity is at « 6 K. However, before a fully 
diamagnetic state could develop, the fluctuating Cooper pair density experiences the pair 
breaking field of the HoNi 5 and hence there is no substantial change in the rising part of 
the M(T) curve. However, at still lower temperature the condensate becomes robust and a 
competition between diamagnetism and ferromagnetism leads to a peak in M(T) at « 3.5 
K. On further lowering the temperature, the diamagnetism of superconducting NbN clearly 
dominates the ferromagnetic response of the system. The appearance of the reentrance in 
trilayers also depends on the thickness {&hn) of each H0M5 layer. As seen in Fig. 1(d), 
when &hn is sufficiently thick (~ 80 nm), the Tq of NbN is greatly suppressed with no 
evidence of reentrance. Also when cLhn is small (~ 20 nm), the Tc of NbN remains robust 
with no sign of entry into the normal state on cooling below T^. In Fig. 1(e) we show 
the variation of Tcurie of trilayer as a function of the total thickness d^ (=<1fi + d^) of 
H0M5 layers. The Tcurie drop by ~ 1.4 K on lowering d^ from 160 nm to 40 nm. From 
these data a very interesting correlation evolves between the Tcurie/Tc of the F-S-F and the 
observation of the reentrance. The later is seen when 0.74 < TcuHe/T^c — 0.92 (see inset of 
Fig. 1(e)). 

In order to get further insight of the reentrant behavior of superconductivity, we have 
measured the resistivity of the sample in the transition region in a magnetic field applied 
along the out-of-plane (Hj_) and in-plane (Hn) directions (perpendicular to bar) of the tri- 
layer. The results of in-plane field measurement are shown in Fig. 2(a). The superconducting 
transition in the H^ configuration is much more sensitive to the field; it becomes signifi- 
cantly broader as we increase the field from to 1400 Oe. An enlarge view of the behavior of 
minimum and peak in R(T) which characterize the reentrant transition is shown in the inset 
of Fig. 2(b) for several values of Hj_. The changes in the minimum have been quantified in 
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terms of AR = R p — R m i n and T min (B.±), and are plotted in the lower inset of Fig. 2(b). The 
AR and T m j„(HjJ drop monotonically with the applied field. The also leads to a large 
number of vortices in the system. The thermally activated flux flow (TAFF) can lead to the 
broadening of the resistive transition in the tail region of the curve. We have analyzed this 
behavior in the framework of TAFF model 34( in the subsequent section. 



The reentrant behavior shows a strong dependence on the current used for the resistivity 
measurements. In Fig. 3(a) we show the R(T) data in zero field for current varying from 1 
yiiA to 2 niA (current density j from 6.6 x 10 2 A/cm 2 to 133.3 x 10 4 A/cm 2 ). The changes in 
the R(T) with current in the temperature regime where the two order parameters compete 
strongly, is quite different from the behavior seen under field (Fig. 2(b)). Here at the 
lowest current density, there is very little evidence of a minimum in the R(T) below the onset 
of superconductivity. We only see a shoulder below which the resistance drops precipitously 
to zero. As the current density is increased from 6.6 x 10 2 A/cm 2 to 6.6 x 10 4 A/cm 2 , 
the minimum in the R(T) becomes pronounced and the T p shows a small shifts to lower 
temperatures (see Fig. 3(b)). A further increase in j to 33.3 x 10 4 A/cm 2 leads to saturation 
of AR. However, for j > 33.3 x 10 A/cm 2 a flattening of the minimum and shift of the Tp 
to lower values is observed (not shown here). The R(T) in the tail portion of the transition 
(T < Tp) becomes progressively broader on increasing the j from 6.6 x 10 2 A/cm 2 . While 
the reasons for this strong current dependence is not understood fully, it is certainly not 
a heating effect. To further clarify the issue of sample heating, we have measured current 
vs. voltage (I-V) curves at different temperatures across the superconducting transition in 
forward and reverse directions. We observe no hysteresis in IVs which could have resulted 
from a thermal lag between the temperature seen by the sensor and the actual temperature 
of the sample. A similar current dependence of the reentrant behavior has been reported by 



Rathnayaka et al. in HoNi 2 B 2 C samples 



271 ] . There authors have reported an interesting 



correlation between the magnitude of the current dependence and the ratio of magnetic and 
superconducting transition temperature. The effect is seen when T m /T c > 0.6. The value of 
T TO /T C in our case is ~ 0.88. Furthermore, our observation of very little reentrant behavior 



27]- 



seen at low currents is also similar to the behavior observed in HoNi 2 B 2 C crystal 

In order to understand how the dissipative behavior of thin NbN in the thermally ac- 
tivated flux flow (TAFF) regime is affected by the formation of F-S-F structures, we have 
compared the R(T) of a plain 10 nm NbN film with that of a F-S and F-S-F structures, hav- 
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ing the same S layer thickness (10 nm). These measurements were performed in Hj_ varying 
in the range of 50 to 1500 Oe. Fig. 4 shows the Arrhenius plots of normalized resistance 
(lnR/R n vs. 1/T) for all three films. We notice a significant broadening of the resistive 
transition with field as we go from the single layer to bilayer and then to trilayer. The 
activation energy for dissipation are calculated using the model used by Hsu and Kapitulnik 



36[ for two dimensional (2D) superconducting film in which the coherence length £(T) and 



penetration depth A(T) are comparable with the film thickness. Since, for NbN £(0) ~ 5 
nm and A(0) ~ 200 nm, and the NbN layer thickness is only ~10 nm, the 2D condition is 
satisfied approximately in this case. Based on the model, the R(T) at the bottom of the 
superconducting transition should follow the relation \ia(R/R N ) = —U (H)/k B T + K(H) 
where Uo(H) is zero temperature activation energy and K(H) is coefficient of linear tem- 
perature correction. Since the low temperature curve is linear, the Uo(H) is obtained by 
calculating the slope of the lnR/R n vs. 1/T curve. The results of this analysis are presented 
in Fig. 5(a). We notice a significant drop in Uo(H) as we go from the plain NbN to F-S 
and then F-S-F. This result can be attributed to the suppression of superconducting order 
parameter at F-S interface due to the proximity coupling with H0M5. Such perturbation on 
the S order parameter is expected to be more in a F-S-F film as both surfaces of the super- 
conducting film see a magnetic layer in this case. This amounts to an effective reduction in 
the thickness of the film by 2£s where £5 is coherence length in S. Similar behavior is also 
observed by Senapati et. al. on high Tc based F-S-F structure where they have shown 
that Uo(H) follow linear relation with thickness of S layer jssj]. 

We further examine the nature of dissipation in all three films by investigating the slope 
of the field dependence activation energy U (H). The inset of Fig. 5(a) shows U (H) in a 
log- log plot. The activation energy of the films at low fields follows the power law relation 
U(H) ~ H- a with a values of « 0.45 , 0.33 and 0.2 for the S, F-S and F-S-F thin films 
respectively. Such variation in the value of the exponent has been observed previously and 
has been explained on the basis of collective vortex creep for a ~ 1, plastic flux creep for 



a ~ 0.5 and single vortex pinning for a ~ 0.1 [37|, l39l44l|. The observed value of a ~ 0.45 
for bare S film suggests the possibility of plastic deformation of flux line lattice in weakly 
pinning vortex liquid. Similar value is also reported by Fogel et. al. at low fields in Mo/Si 



multilayer system 



421 ] . The presence of HoNis layer appears to enhance vortex pinning, 



leading to a single vortex depinning dynamics. 
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We now address the strong measuring current dependence of the dissipation in the tail 



region of the superconducting transition in the F-S-F film [43]. The data in the tail region 
of the transition are plotted as a function of 1/T in the insert of Fig. 5(b). The behavior of 
activation energy for a fixed j value is shown in the main panel. The Uq(H) remains nearly 
constant for j < 10 x 10 4 A/cm 2 but beyond this value a sharp drop in the activation energy 
is seen. This suggest a transition to the collective creep regime of vortex motion at high 
current. 

We close the discussion by proposing a possible scenario for the observation of reentrance 
in superconductivity in our magnetic thin film heterostructures. It is important to point 
out here that our observation of reentrant superconductivity as a function of temperature 
appears different from the so called reentrant behavior recently reported in Nb-Cu4iNi 5 g 



bilayers as a function of CuNi layer thickness |6|,|7j]. Such observations have been made earlier 
as well {4, 5|. This arises because the superconducting pair amplitude (PA) is oscillating 
inside the F, and depending on the thickness of the F layer [j], Q], the PA at the interface 
region of the F can be very large or nearly zero, leading to a higher or lower T c respectively. A 
sudden entry into (or out of) the normal state can also occur if the magnetization vectors of 
the F layers transit from a parallel to antiparallel configuration as a function of temperature. 



These spin accumulation 



44| and Cooper pair averaging effect have been addressed in the 



past [45] . However, such magnetization reorientation transition is more likely to occur in a 
magnetic field than by temperature. Recently Wu, Vails and Halterman [25] have proposed 
an interesting possibility of observing a temperature dependent reentrant superconductivity 
in a F-S film if the F layer has a spiral magnetic order, such as that exists in Holmium. The 
magnetic structure of HoNi 5 is not fully understood, although magnetization measurements 
show strong anisotropy with a-b plane of hexagonal unit cell being the easy plane. The H0M5 
films grown on MgO substrate shows inhomogeneities in magnetization due to different 
orientation of nanostructures, as we have discussed in previous report {m]. It is possible 
that the inhomogeneous magnetization of H0M5 allows a reentrant behavior in the Wu, Vails 



and Halterman sense 



25] . This effect is presumably accentuated by the presence of HoNi 5 



on both sides of the thin NbN layer. 

A more intuitive explanation for the RES can, however, be given in terms of proxim- 
ity effect in a normal (N) and ferromagnetic metal. The decay of superconducting order 

1 /2 

parameter in N and F layers is given by the length scale — {hD N /2Trk B T C s) and 
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C F = (hD F /27rK B T cs ) 1/2 , respectively 151 

. Here £jv an d £f are coherence length, and D N 
and Dp are diffusion coefficients in the normal and ferromagnetic metal respectively. The 

1 /2 

other relevant length scale in the problem is £p = (hDp/E ex ) 1 , which defines the decay of 
the S order parameter in the F layer. It is measure of the length scale over which supercon- 
ductivity is induced in the F layer. The exchange splitting of the conduction band E ex is 
related to the exchange integral I and magnetic moment fip as E ex = I\xp [46j. Since E ex 
goes to zero at Tcurie, £f diverges on warming the sample toward the magnetic ordering 
temperature. This divergence is however cut off by £at. This is the unique feature of our 
experiment because Tcurie < Tc. As long as Tcurie < T < Tc, the superconductivity in 
the effective thickness of ds + 2£at where ds is the thickness of superconducting layer, leads 
to a drop of resistance. But as soon as T ~ Tc urie this effective thickness starts decreasing 
not only because £p < £jv but also because of the pair breaking effects of the magnetic layer 
inside the superconducting film. The consequence of this would be a rise in the resistance. 
We have estimated the diffusion coefficient Dp of H0M5 from the specific heat data 47J, |48] 



and measured resistivity of 40 nm thick HoNi 5 film, which is ~ 280 /ificm. Since the elastic 
mean free path for this value of resistivity is already of the order of interatomic distance, 
we expect a marginal increase of resistivity of the thinner H0N15 film due to size effect 49]. 
Band structure calculations of E ex for HoNi 5 5(3] suggest a value of 0.5 eV, which yields 
£f(0) ~ 0.3 nm. Since for a strong ferromagnet like Fe, £p is only 1 nm [51], we believe 
that this calculation overestimates E ex . The £n on the other hand is 2.7 nm. As the corre- 
lation length in HoNi 5 reduces on going through the magnetic transition, a reentrance into 
the normal state is expected. However, at still lower temperatures, the order parameter in 
the S film becomes robust enough to short circuit the H0M5 films and the resistance would 
start approaching zero value. These semi quantitative arguments are consistent with the 
results in Fig. 1. 



CONCLUSION 



In summary, we have carried out a detail study of electronic transport in F-S and F-S-F 
thin film where the F and S order parameters are of comparable strength. In the F-S-F 
trilayer we see a robust reentrant superconductivity over a critical range of the ratio of mag- 
netic and superconducting transition temperature of the F-S-F trilayer (0.74 < Tc-wie/Tc 
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< 0.92). The RES seen here is different from the more common reentrant behavior observed 
as a function of sample geometry. Many features of the RES reported here are similar to 
those found in bulk magnetic superconductors consisting of ternary and quaternary alloys 
of rare earths (29), 3p| . 



ACKNOWLEDGMENTS 

This research has been supported by a grant from Department of Information Technology 
(DIT), India. Gyanendra Singh acknowledges financial support from the Council of Scientific 
and Industrial Research (CSIR), Government of India. RCB acknowledges financial support 
from J.C. Bose Fellowship of the Department of Science and Technology, India. 



Electronic address: rcb@iitk.ac.in, rcb@nplindia.org 



[1] A. I. Buzdin, Rev. Mod. Phys. 77, 935 (2005). 

[2] I. F. Lyuksyutov and V. L. Pokrovsky, Adv. Phys. 54, 67 (2005). 

[3] E. A. Dernier, G. B. Arnold, and M. R. Beasley, Phys. Rev. B 55, 15174 (1997). 

[4] Th. Muhge, K. Theis-Brohl, K. Westerholt, H. Zabel, N. N. Garifyanov, Yu. V. Goryunov, I. 

A. Garifullin, and G. G. Khaliullin, Phys. Rev. B 57, 5071 (1998). 
[5] C. Cirillo, S. L. Prischepa, M. Salvato, C. Attanasio, M. Hesselberth, and J. Aarts, Phys. Rev. 

B 72, 144511 (2005). 

[6] V. Zdravkov, A. Sidorenko, G. Obermeier, S. Gsell, M. Schreck, C. Muller, S. Horn, R. Tidecks, 
and L. R. Tagirov Phys. Rev. Lett. 97, 057004 (2006). 

[7] V. I. Zdravkov, J. Kehrle, G. Obermeier, S. Gsell, M. Schreck, C. Muller, H.-A. Krug von 
Nidda, J. Lindner, J. Moosburger- Will, E. Nold, R. Morari, V. V. Ryazanov, A. S. Sidorenko, 
S. Horn, R. Tidecks, and L. R. Tagirov Phys. Rev. B 82, 054517 (2010). 

[8] I. A. Garifullin, D. A. Tikhonov, N. N. Garifyanov, L. Lazar, Yu. V. Goryunov, S. Ya. Khleb- 
nikov, L. R. Tagirov, K. Westerholt, and H. Zabel, Phys. Rev. B 66, 020505(R) (2002). 

[9] I. A. Garifullin, D. A. Tikhonov, N. N. Garifyanov, M. Z. Fattakhov, L. R. Tagirov, K. Theis- 
Brohl, K. Westerholt, and H. Zabel, Phys. Rev. B 70, 054505 (2004). 
[10] Y. Obi, M. Ikebe, and H. Fujishiro, Phys. Rev. Lett. 94, 057008 (2005). 



10 



[11] G. Verbanck, C. D. Potter, V. Metlushko, R. Schad, V. V. Moshchalkov, and Y. Bruynseraede, 

Phys. Rev. B 57, 6029 (1998). 
[12] J. S. Jiang, D. Davidovic, Daniel H. Reich, and C. L. Chien, Phys. Rev. Lett. 74, 314 (1995). 
[13] Th. Muhge, N. N. Garifyanov, Yu. V. Goryunov, G. G. Khaliullin, L. R. Tagirov, K. Wester- 
holt, I. A. Garifullin, and H. Zabel, Phys. Rev. Lett. 77, 1857 (1996). 
[14] L. V. Mercaldo, C. Attanasio, C. Coccorese, L. Maritato, S. L. Prischepa, and M. Salvato, 

Phys. Rev. B 53, 14040 (1996). 
[15] T. Kontos, M. Aprili, J. Lesueur, and X. Grison, Phys. Rev. Lett. 86, 304 (2001). 
[16] V. V. Ryazanov, V. A. Oboznov, A. Yu. Rusanov, A. V. Veretennikov, A. A. Golubov, and J. 

Aarts, Phys. Rev. Lett. 86, 2427 (2001). 
[17] V. A. Oboznov, V. V. Bolginov, A. K. Feofanov, V. V. Ryazanov, and A. I. Buzdin, Phys. 

Rev. Lett. 96, 197003 (2006). 
[18] J. W. A. Robinson, S. Piano, G. Burnell, C. Bell, and M. G. Blamire, Phys. Rev. Lett. 97, 

177003 (2006); Phys. Rev. B 76, 094522 (2007). 
[19] T. S. Khaire, M. A. Khasawneh, W. P. Pratt Jr., and N. O. Birge, Phys. Rev. Lett. 104, 

137002 (2010). 

[20] J. W. A. Robinson, J. D. S. Witt, and M. G. Blamire, Science 329, 59 (2010). 

[21] F. S. Bergeret, A. F. Volkov, and K. B. Efetov, Rev. Mod. Phys. 77, 1321 (2005). 

[22] P. V. Leksin, N. N. Garifyanov, I. A. Garifullin, Ya. V. Fominov, J. Schumann, Y. Krupskaya, 

V. Kataev, O. G. Schmidt, and B. Buchner, Phys. Rev. Lett. 109, 057005 (2012) 
[23] J. Zhu, I. N. Krivorotov, K. Halterman, and O. T. Vails, Phys. Rev. Lett. 105, 207002 (2010). 
[24] Gyanendra Singh, P. C. Joshi and R. C. Budhani, Physica C 472, 44 (2012). 
[25] Chien-Te Wu, Oriol T. Vails and Klaus Halterman, Phys. Rev. Lett. 108, 117005 (2012). 
[26] H. Eisaki, H. Takagi, R. J. Cava, B. Batlogg, J. J. Krajewski, W. F. Peck, Jr., K. Mizuhashi, 

J. O. Lee, and S. Uchida, Phys. Rev. B 50, 647 (1994). 
[27] K. D. D. Rathnayaka, D. G. Naugle, B. K. Cho, and P. C. Canfield, Phys. Rev. B 53, 5688 

(1996). 

[28] L. C. Gupta, Advances In Physics 55, 691 (2006). 

[29] W. A. Fertig, D. C. Johnston, L. E. DeLong, R. W. McCallum, M. B. Maple, and B. T. 

Matthias, Phys. Rev. Lett. 38, 987 (1977). 
[30] G. W. Crabtree, F. Behroozi, S. A. Campbell, and D. G. Hinks, Phys. Rev. Lett. 49, 1342 

11 



(1982). 

[31] Gyanendra Singh, P. C. Joshi and R. C. Budhani, J. Appl. Phys. 109, 113915 (2011). 
[32] K. Senapati, N. K. Pandey, Rupali Nagar, and R. C. Budhani, Phys. Rev. B 74, 104514 
(2006). 

[33] D. Gignoux, A. Nait-Saada, A. Perrier de la Bthie, J. Phys. Coll. 40 (C5) (1979) 188. 
[34] P. W. Anderson and Y. B. Kim, Rev. Mod. Phys. 36, 39 (1964). 

[35] P. H. Kes, J. Aarts, J. van der Beek, and J. A. Mydosh, Supercond. Sci. Technol. 1, 242 
(1989). 

[36] J. W. P. Hsu and A. Kapitulnik, Phys. Rev. B 45, 4819 (1992). 

[37] K. Senapati and R. C. Budhani, Phys. Rev. B 70, 174506 (2004). 

[38] W. R. White, A. Kapitulnik, and M. R. Beasley, Phys. Rev. Lett. 70, 670 (1993). 

[39] T. T. M. Palstra, B. Batlogg, L. F. Schneemeyer, and J. V. Waszczak, Phys. Rev. Lett. 61, 

1662 (1988).; T. T. Palstra, B. Batlogg, R. B. van Dover, L. F. Schneemeyer, and J. V. 

Waszczak, Phys. Rev. B 41, 6621 (1990). 
[40] J. T. Kucera, T. P. Orlando, G. Virshup, and J. N. Eckstein, Phys. Rev. B 46, 11004 (1992). 
[41] O. Brunner, L. Antognazza, J.-M. Triscone, L. Miville, and O. Fischer, Phys. Rev. Lett. 67, 

1354 (1991). 

[42] N. Ya. Fogel, V. G. Cherkasova, O. A. Koretzkaya, and A. S. Sidorenko, Phys. Rev. B 55, 85 
(1997). 

[43] Y. Abulafia, A. Shaulov, Y. Wolfus, R. Prozorov, L. Burlachkov, Y. Yeshurun, D. Majer, E. 

Zeldov, H. Wuhl, V. B. Geshkenbein, and V. M. Vinokur, Phys. Rev. Lett. 77, 1596 (1996). 
[44] P. G. de Gennes, Rev. Mod. Phys. 36, 225 (1964). 
[45] Y. Fominov thesis, |arXiv:cond-mat / 031 1359 , 

[46] J. Aarts, J. M. E. Geers, E. Bruck, A. A. Golubov, and R. Coehoorn, Phys. Rev. B 56, 2779 
(1997). 

[47] S. G. Sankar, D. A. Keller, R. S. Craig, W. E. Wallace, And V. U. S. RAO, J. Solid State 

Chemistry 9,78 (1974). 
[48] C. A. Bechman, W. E. Wallace and R. S. Craig J. Phs. Chem. Solids. 3 463 (1974). 
[49] P. A. Lee and T. V. Ramakrishnan, Rev. Mod. Phys. 57, 287 (1985). 
[50] S. K. Malik, F. J. Arlinghaus and W. E. Wallace, Phys. Rev. B 25, 6488 (1982). 
[51] P. V. Leksin, N. N. Garifyanov, I. A. Garifullin, J. Schumann, V. Kataev, O. G. Schmidt, and 

12 



B. Buchner, Phys. Rev. Lett. 106, 067005 (2011). 



13 




FIG. 1: (Color online) (a) X-ray reflectivity of HoNis/NbN bilayers shown as open circles with 

NbN thickness of 10 and 30 nm. Solid line is the best fitting of genetic algorithm. The upper 

inset shows the geometry of patterned sample for four probe measurements. The lower inset shows 

resistance vs. temperature of 10 nm thick NbN single layer (B) and H0M5 /NbN/HoNis trilayers 

with NbN thickness 10 nm (A) and 30 nm (C). (b) Larger view of R(T) for trilayer with 10 nm 

NbN. The thickness of each H0M5 layer in (b) and (c) is 50 nm. A colored band between 5.2 to 

5.46 K indicates the region of reentrant superconductivity, (c) Magnetization vs. temperature of 

trilayer with 10 nm thick NbN, measured in 100 Oe out-of-plane field after zero field cooling of the 

film. The inset shows M vs. T between 3 to 7 K to highlight the magnetic ordering temperature of 

the film. It also shows M vs T for single layer H0M5 thin film measured after ZFC and applying 

the 100 Oe field in out-of-plane direction. A small but distinct ferromagnetic contribution to M(T) 

is seen at T > Tcurie which we believe comes from uncorrelated 4f spins of Ho. (d) R(T) of 

HoNi5/NbN/HoNi5 trilayers at zero field with five different thicknesses of H0M5 ~20, 40, 50, 60, 

80 nm. The NbN thickness was kept constant to ~10 nm for all the samples, (e) Variation of the 

magnetic ordering temperature of the trilayer as a function of the total thickness of HoNis layer. 

Inset show the dependence of Tcurie/^c on d^. The shaded area range where reentrance is seen. 
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FIG. 2: (Color online) (a) Temperature dependent resistance at constant in-plane magnetic field 
varying from to 1400 Oe. The inset shows larger view of the data for reentrant region, (b) R(T) 
with field applied along out-of-plane direction. A 100 /iA current is applied along the bar for these 
measurements. The upper inset shows the clear view of reentrant behavior as a function of H^. 
The peak and dip in R(T) is marked by Hp and R m j n respectively. The lower inset shows the 
calculated value of AR (= R p - R m j n ) and T m j n as a function of Hj_ 
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FIG. 3: (Color online) (a) Temperature dependent resistance measured in zero magnetic field with 
constant current density j between 6.6 x 10 2 A/cm 2 to 33.3 x 10 4 A/cm 2 (current 1 /iA to 500 //A). 
The upper inset shows magnified view of R(T) where normal state reentrant behavior is observed. 
The behavior of critical temperature (Tc) with applied current is plotted in the lower inset, (b) 
Peak position (Tp) and AR (= R p - R m in) as a function of current density calculated from panel 
(a). 
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FIG. 4: (Color online) Normalized resistance vs. 1/T (Arrhenius plots) for (a) NbN(10 nm)/MgO 
(b) NbN(10 nm)/HoNi 5 (50 nm)/MgO bilayer and (c) HoNi 5 (50 nm)/NbN(10 nm)/HoNi 5 (50 
nm)/MgO trilayer at various value of a magnetic field applied in the out-of-plane geometry (Hj_). 
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FIG. 5: (Color online) (a) Activation energy as a function of magnetic field applied perpendicular to 
plane of NbN(10 nm) /MgO, NbN(10 nm) /HoNi 5 (50 nm) /MgO bilayer and HoNi 5 (50 nm) /NbN(10 
nm)/HoNi5 (50 nm)/MgO trilayer. The inset shows plot in log-log scale. Solid line is the best fit 
to U{H) ~ H~ a . (b) Current dependence effective activation energy of the trilayer extracted from 
the Arrhenius plots shown in the inset of the figure. The value of 5 is calculated from the best 
fitting of U e // a j _<5 . 
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